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We study the mechanism of the triplet superconductivity (TSC) in Sr2Ru04 based on the multiorbital Hub¬ 
bard model. The electronic states are studied using the recently developed renormalization group method 
combined with the constrained random-phase-approximation, called the RG-l-cRPA method. Thanks to the 
vertex correction (VC) for the susceptibility, which is dropped in the mean-field-level approximations, strong 
orbital and spin fluctuations at Q « (Stt/S, 27r/3) emerge in the quasi one-dimensional Fermi surfaces (FSs) 
composed of -|- dy^ orbitals. Due to the cooperation of both fluctuations, we obtain the triplet supercon¬ 
ductivity in the Eu representation, in which the superconducting gap is given by the linear combination of 
{Ax{k), Ay{k)) ~ (sin Sfcj;, sin Sfcj,). Very similar results are obtained by applying the diagrammatic cal¬ 
culation called the self-consistent VC method. Thus, the idea of “orbital-l-spin fluctuation mediated TSC” is 
confirmed by both RC-l-cRPA method and the self-consistent VC method. We also reveal that a substantial 
superconducting gap on the d^^y-orbital FS is induced from the gaps on the quasi one-dimensional FSs, in con¬ 
sequence of the large orbital-mixture due to the 4d spin-orbit interaction. 

PACS numbers: 74.20.Rp, 71.27.-ta, 74.25.Jb, 74.70.Pq 


I. INTRODUCTION 

Sr2Ru04 is an unconventional smerconductor with the 
transition temperature Tc = 1 . 5 K This material has 

been attracting great attention since the spin triplet supercon¬ 
ductivity (TSC) is indicated by the NMR study 01 . The chiral 
p-wave {px + ipy) TSC, which is analogous of the A-phase of 
the superfluid ^He, had been predicted in Ref. ist]. However, 
in contrast to superfluid ^He, the paramagnon mechanism is 
not realized in Sr2Ru04 since no ferro-magnetic fluctuations 
are observed. Instead, strong antiferro-magnetic (AFM) fluc¬ 
tuations with Q Ri ( 2 W 3 , 27 r/ 3 ) are observed by neutron 
scattering spectroscopy 0 . Since the AFM fluctuations usu¬ 
ally drive the spin singlet superconductivity (SSC), the mech¬ 
anism of the TSC in Sr2Ru04 has been a long-standing prob¬ 
lem in strongly correlated electron systems. 

Figures [T](a) and (b) show the (dxz, dj,z)-orbital bands and 
the corresponding quasi-one-dimensional (qlD) Fermi sur¬ 
faces (FSs), FSa and FS/ 3 . In Sr2Ru04, the nesting of the 
qlD FSs gives the experimentally observed AFM fluctuations 
at Q r: ( 27 r/ 3 , 27 r/ 3 ). In addition, there is a da^y-orbital band 
that compose the two-dimensional ( 2 D) FS, FS7. (see Figs. 
| 4 ](a) and (b).) If the spin-orbit interaction (SOI) is neglected, 
the (a, ,8)-bands and 7-band are coupled only via the Coulomb 
interaction. Therefore, the superconductivity is expected to be 
realized mainly in either the qlD bands (|Aq, p| 3> |A-y|) or 
the 2 D band (|Aq,^p| <?C |A..y|). 

The mechanisms of the TSC originating mainly from the 
2 D band had been proposed in Refs. ItI- H^ : Nomura and 
Yamada explained the TSC state by using the higher-order 
perturbation theory ||3l, which is the natural development of 
the Kohn-Luttinger mechanism d. Recently, a three-orbital 
Hubbard model had been studied using a 2 D renormalization 
group (RG) method |[^, and it was claimed that the p-wave 
gap is realized on the FS7 accompanied by the development 
of spin fluctuations at q = (O.IQtt, O.IQtt). Also, charge- 


fluctuation-mediated TSC was discussed by introducing the 
inter-site Coulomb interaction dst]. 

On the other hand, it is natural to expect that the TSC is 
closely related to the strong AFM fluctuations in the qlD 
FSs at q ~ Q. The TSC originating from the qlD FSs had 
been discussed by applying the perturbation theory dQ 
and random-phase-approximation (RPA) dQ. Takimoto 
discussed the orbital-fluctuation-mediated TSC using the RPA 
under the condition U' > U, where U (If) is the intra-orbital 
(inter-orbital) Coulomb interaction 01 ^ . However, in the 
RPA, the SSC is obtained under the realistic condition U > U' 
due to strong AFM fluctuations. When the spin fluctuation is 
Ising-like, the TSC may be favored since the pairing inter¬ 
action for the SSC is reduced 111 ^. T he charge-fluctuation- 
mediated TSC was also discussed fl^ . However, these theo¬ 
ries did not clearly explain why the TSC overwhelms the SSC 
in Sr2Ru04 despite the presence of strong AFM fluctuations. 

To And out the origin of the TSC in Sr2Ru04, many exper¬ 
imental efforts have been devoted to determine the gap struc¬ 
ture, such as the tunnel junction ifp^ . ARPES, and quasiparti¬ 
cle interference measurements. Recently, large superconduct¬ 
ing gap with 2 | A| r 5 Tc was observed by the scanning tun¬ 
neling microscopy measurements iH. The observed large 
gap would be that on the qlD FSs, since the tunneling will 
be dominated by the (d^^z, dyz)-orbitals that stand along the 
z-axis, as clarified in the double-layer compound Sr3Ru207 
Therefore, it is an important challenge to establish the 
theory of the TSC based on the qlD-band Hubbard model, by 
applying an advanced theoretical method. 

In this paper, we study the multiorbital Hubbard models for 
Sr2Ru04 with realistic parameters {U > U'). To analyze 
the many-body effect beyond the mean-fleld-level approxi¬ 
mations, we apply both the RG-tcRPA method developed in 
Ref. B in Sec.lfl and the SC-VCs method in Refs. mf. B 
in Sec.lnl] Thanks to the vertex correction (VC) for the sus¬ 
ceptibility, which is dropped in the RPA, strong orbital and 
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FIG. 1: (color online) (a) Band structure and (b) FSs of the two- 
orbital model. Q « {2n/3,2n/3) is the nesting vector, (c) X’Iq) 
and (d) x^{q) of Ihs qlD-band model obtained by the RG-l-cRPA 
method (Aq = 1) for U = 3.5, J/U = 0.035 and T = 0.02. 


spin fluctuations at Q Ri (27r/3, 27r/3) emerge in the qlD 
bands, and it is confirmed that the TSC is efficiently induced 
by the cooperation of both fluctuations for wide parameter re¬ 
gion. The idea of the “spin-Horbital fluctuation mediated TSC 
state” will be realized not only in Sr 2 Ru 04 , but also in other 
triplet superconductors. In Sec. IIIIBI we discuss the effect 
of the SOI on the superconductivity, and reveal that the large 
gap on the FS 7 is induced from the FS/3 due to the SOI- 
induced “orbital-mixture” between FS/3 and FS 7 . This effect 
will be important to understand the experimental approximate 
T -linear behavior of C/T below Tc in Sr 2 Ru 04 . 


II. RENORMALIZATION GROUP + CRPA METHOD 

In this section, we analyze this model by applying the RG 
combined with the constrained RPA (RGh-cRPA) This 
method is very powerful to calculate the higher-order many- 
body effects systematically and in an unbiased way. In the 
RGh-cRPA method, we divide the lower-energy region (| E\ < 
An ) of the Brillouin zone into N patches as done in Refs. 11231 - 
1^ and perform the RG analysis. The contributions from the 
higher-energy region (|i3| > Aq) are calculated by the cRPA 
method with high numerical accuracy, and incorporated into 
the initial vertex functions 1^ . (The conventional patch-RG 
method ||23| - I^ is recovered when Aq > Wband-) Although 
the initial vertex functions are very small, they play decisive 
roles for the fixed point of the RG flow. 

A. Analysis of the Two-Orbital Model 

First, we study the two-orbital Hubbard model, which de¬ 
scribes the quasi-lD FSs of Sr 2 Ru 04 . The kinetic term is 
given by 


1,2 

( 1 ) 

k,(T l,m 

where the orbital indices l,m = 1 and 2 refer to dxz- and 
dy 2 -orbitals, respectively. In the present model, = 


— 2tCOSkx — 2tnn cosky, = —2t cosky — 2tanCOSkx, 

and = At' sinkx sinky. Hereafter, we set (t,tnn,t') = 
(1, 0.1, 0.1), and fix the Ailing as n = 4 • (2/3) = 2.67, which 
corresponds to the Ailing of the qlD FSs of Sr 2 Ru 04 . We 
also introduce the on-site Coulomb interactions U, U', and 
put the exchange and Hund’s couplings J = J' = {U — U')/2 
throughout the paper. 

Here, we analyze this model by applying the RG-rcRPA 
method j^ . We use N — 6A (32 patches for each FS) in 
the present study, and it is verified that the results of N = 128 
are almost unchanged. First, we calculate the susceptibilities 
using the RG-rcRPA: The charge (spin) susceptibility is given 
by 




,( 2 ) 

(3) 


where q = and l,l',m,m' are d orbitals. The 

quadrupole susceptibility with respect to Oa; 2 _j ^2 = Uxz—nyz 
and the total spin susceptibility are respectively given as 


X^{q) = ( 4 ) 

l,m 

X"{q) = '^Xll-m,miQ)- ( 5 ) 

l,m 

Figures [T](c) and (d) show the obtained x^{q) x^(9)> re¬ 

spectively, by the RG-rcRPA method (Aq = 1) for (7 = 3.5 
and J/U = 0.035 at T = 0.02. Both susceptibilities have 
the peak at Q ss (27r/3, 27r/3), which is the nesting vector of 
the present FSs. The shape of x^{q) essentially equivalent 
to that of the RPA, by putting U = 2.2 and J/U = 0.035. 
However, x^io) in the RPA is quite small when J > 0 
ll^ 1^ . Therefore, the enhancement of x^{q) in Fig- [11(d) 
originates from the many-body effect beyond the RPA. The 
natural candidate is the Aslamazov-Larkin (AL) type VC for 
x'^(q), Ar°(g), whose analytic expression is given in Ref. li^ . 
Since X‘'{q) ^ Aal( 9; fc)^x®(^)x®(fc + O') for sim¬ 

plicity, X^{q) takes large value at q = 0 and 2Q when x*(fc) 
is large at k = Q. AAhiq; k) is the three-point vertex com¬ 
posed of three Green functions J^. In the present model, 
2Q Ri <5 in the first Brillouin zone. Thus, with the aid of the 
VC and the nesting of the FSs, the enhancement of X^{Q) in 
Fig. [11(d) is realized. 

Figure|2](a) shows x^(Q) nnd x'^i.Q) ns functions of J/U 
at T = 0.02, obtained by the RGh-cRPA method with Aq = 1. 
For each value of U, x^(Q) (x^{Q)) decreases (increases) 
with J/U, and they are equal at {J/U)c ~ 0.035. We stress 
that {J/U)c is negative in the RPA since the VC is totally 
dropped. In the case of Aq = 3 shown in Fig.|2(b), the value 
of {J/U)c increases to ^ 0.09 at (7 ~ 5.5, indicating that 
importance of the VC due to higher energy region. To check 
this expectation, we include the constrained AL term (cVC) 
in addition to the cRPA 1I20I] . The obtained results are shown 
in Fig. |2](c). It is verified that {J/U)c increases to 0.08 at 
U = 3.6. (x' 5(Q) in Fig. id (c) is approximately given by 
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FIG. 2: (color online) X*(Q) X^{Q) functions of J/U given 

by the RG+cRPA method for (a) Aq = 1 (i7 = 3.0 ~ 3.8) and (b) 
Ao = 3 (U — 4.0 ~ 5.5). (c) x^iQ) x^iQ) for = 1, 
by including the cVC. (d) Obtained phase diagram for Aq = 1 and 
Ao = 3 without cVC. 


shifting {Q) in Fig-ll](a) horizontally by +0.02 ^ +0.05.) 
The values of {J/U)c obtained by Figs.|2](a)-(c) are summa¬ 
rized in Fig. |2](d). We find that {J/U)c ~ 0.15 in the self- 
consistent VC method with the self-energy correction (SC- 
VCe method) HI. 

Although the value of {J/U)c is underestimated at Aq = 1, 
the obtained x^{q) x'^{q) at Aq = 1 are reliable, since 

the higher-energy processes can be calculated with high nu¬ 
merical accuracy 12^ . Hereafter, we perform the RGh-cRPA 
method with Aq = 1, by using smaller J/U (~ 0.04) to 
compensate for the absence of the higher-energy VCs. Fig¬ 
ure |3] (a) shows the T -dependences of X®(Q) and x^{Q) 
given by the RGh-cRPA method (Aq = 1) for U = 3.8 and 
J/U = 0.04; Both of them are strongly renormalized from 
the RPA results. In the RPA, Xrpa(Q) diverges at T sa 0.4, 
at which Xrpa(Q) remains very small. In highly contrast, in 
the RGh-cRPA method, the relation x*(Q) ~ x'^i.Q) holds 
for wide temperature range. 

We also calculate the TSC and SSC susceptibilities using 
the RGh-cRPA method; 

xl^s) = ljyr{Bl^s)i^)B,^sm). ( 6 ) 

where ^ = «’/3is 

the band index, and (q) is the odd (even) parity gap func¬ 
tion. The obtained xf(^) is shown in Fig.|3(a), by optimizing 
the functional form of (q) numerically ll2^ . Since xf(^) 

diverges at T = Tc, the strong development of Xt ^ at T ss 
0.02 means that the TSC is realized. This TSC state belongs 
to the two-dimensional i?u-representation, (A^(q), A^(q)). 
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FIG. 3; (color online) (a) T-dependences of X*(Q)> X^{Q)^ Xs*^ 
and Xt^ for U = 3.8 and J/U — 0.04 (Aq = 1). (b) gap 

functions on FS^, A^{9) (/r = a,/3) obtained by the RG. The 
relation Af oc sin3fca; holds approximately. N = 128 patches 
are used, ( c) The magnitud e of the chiral (or helical) gap state 
I A''I = . (d) Schematic explanation for the 

sin3fcx-type TSC due to orbitalH-spin fluctuations at q = Q. Solid 
lines (broken lines) are the necessary (accidental) nodes. The posi¬ 
tions of nodes (A^ = 0) in (b) are shown by crosses, (e) Aig and (f) 
Big SSC gap functions. 


The obtained A^ on the FSs when Xt^ 60 are shown 
in Fig. [3] (b), where 9 is the angle of the Fermi momentum 
shown in Fig. [T](b). The necessary nodes = 0 are on 

the lines qx{y) = 0, ±7r. Very similar TSC gap is obtained for 
J/U < 0.08 by taking the cVC into account with Ag = 1. Be¬ 
low Tc, the BCS theory tells that the chiral or helical gap state 
with the gap amplitude |A^| = y/(A^)^ + [AyY, which is 
shown in Fig. |3 (c), is realized to gain the condensation en¬ 
ergy. 
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To understand why the TSC state is obtained, it is useful to 
analyze the linearized gap equation; 


AfA;^(g) = 


dq' 


X ln( 1 . 13 a;c/r), 


{q.q’)K (?') 


(7) 


where a = t or s. Af' is the eigenvalue, {q, q') is the 
pairing interaction, and ujc is the cut-off energy of the interac¬ 
tion. As shown in Fig.[T](b), the inter-band interaction (/r = a, 
fi' = /?) with q~ q' = Q is approximately given by the intra¬ 
orbital interaction given as 




^a^|Af(q;q')Px?(g-g') 

+ Ca^|A^(q;q')l^xf(9-g'): (8) 



where (6t,ct) = (-1,-1) and (6s, Cs) = (3,-1), and 
xt’‘'iQ) = XiS-i,i{Q)- (Note that x?(Q) ~ x%Q)/^ 
and xliQ) ~ X^{Q)/y since xtiQ) > Xf 1,2^2 (Q) and 
XiiQ) ~ “Xi,i; 2.2(Q) near the critical point 11^ .') is 
the VC for the gap equation, which we call A-VC in Ref. 
d^. The AL-type diagram for the charge channel is given by 
AF(q;g') ~ l+TJ 2 k^AL{q-q']k)Gik)x"ik+q)x"{k-q'), 
which is strongly enlarged for q — q' k, Q, and the orbital- 
fluctuation-mediated pairing is favored 0^ . The merit of 

the RGh-cRPA method is that the the AL-type A-VC is auto¬ 
matically produced in calculating the pairing susceptibility in 
Eq. (|6]l. 

In the RPA with J > 0, the TSC cannot be achieved be¬ 
cause of the relation XiiQ) ^ XiiQ) A'^’® = 1 in 
the RPA: In this case, spin-fluctuation-mediated SSC is ob¬ 
tained since |T4*| Ri (C/^/2){3|Afpxf ~ |AfPxF} N three 
times larger than |V/| Ri (C/V2){|Af-f |Afpxf}. In 
the present RG-rcRPA method, in contrast, the relationship 
xf(Q) ~ xF {Q) is realized, and therefore the triplet interac¬ 
tion \V^ I can be larger than |. We verified that the TSC sus¬ 
ceptibility overcomes the s-wave and d-wave SSC susceptibil¬ 
ities when (Q) is comparable or larger than X* {Q), which 
is realized in wide parameter range in the present RG-rcRPA 
method as shown in Figs.|2](a)-(c). 

Using Fig. |3](d), we explain the gap structure of the TSC 
state induced by orbital-rspin fluctuations si q ~ Q. In ad¬ 
dition to the necessary nodes shown by solid lines, accidental 
nodal lines appear around ~ ±7r/3 and kx Rs ±27r/3: The 
reason is that A“(q) and A^{q') tend to have the same sign 
for q — q'^iQ due to large attractive interaction by V^{q,q'). 
For this reason, the relation A^(q) ~ sinSkx in Fig.|3](b) is 
satisfied in the Eu-type TSC state. 

In Fig. [3 (a), also develops at low temperatures; Fig¬ 
ures |3](e) and (f) show the obtained Aig and Big SSC gap 
structures, which give the first and the second largest 
Both SSC states with sign reversal are mainly caused by spin 
fluctuations, and Aig state is slightly stabilized by the or¬ 
bital fluctuations. The Aig state in Fig. |3](e) dominates the 
TSC state when x^(Q) ^ X^iQ)^ which is realized for 
J/U > 0.05 in Fig. 12(a). 


FIG. 4: (color online) (a) Band structure, (b) three FSs, and (c) 
DOS in the three-orbital model for Sr 2 Ru 04 without the SOI. In the 
present model, ALji(0)/Atotai(0) = 0.42, which is consistent with 
the value of the band calculation (0.47). 


B. Analysis of the Three-Orbital Model 

In the previous subsection, we showed that the 
“spin-Horbital fluctuation mediated TSC state” emerges 
by analyzing the two-orbital Hubbard model. In order to ver¬ 
ify this result, we study a more realistic three orbital Hubbard 
model for Sr2Ru04: We add the 7-band given by the dxy- 
orbital, = — 2 t^{cos kx + cos ky) — dfg cos kx cos ky + E 3 , 
to the present two-orbital model HU]- We put tg = 0.86, 
fg = 0.36, E 3 = O.OI, Hxz + riyz = 2.67 and Uxy = 1.28. 
In Fig. m we show the (a) band structure, (b) FSs, and (c) 
density-of-states (DOS) in the present three-orbital model in 
the absence of the SOI. 

We analyze this three-orbital Hubbard model by applying 




FIG. 5: (color online) (a) x^ (q) and (b) x^ (q) given by applying the 
RG-l-cRPA method to the three-orhital model, (c) The TSC gap func¬ 
tions as function of 6. (d) Obtained temperature dependence of the 
TSC (SSC) susceptibility, Xt(2), together with x*^(Q) and x‘iQ)- 
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the RG+cRPA method. Here, we take the renormalization of 
the Green function due to (S = a,/ 3 , 7 ) according to 
Ref. S, where Zs is the mass-enhancement constant of 6 - 
band, and we put Z^/Z^^p = 1.4 (Z„ a = 1) as observed by 
ARPES and dHvA measurements 11^ . and as theoretically 
obtained by the SC-VCs method (see Fig. | 6 ](c)). Figures |5] 
(a) and (b) show the obtained x^{q) respectively. 

We put Aq = 1, U = 3.9 and J/U = 0.02. Thus, the strong 
orbital and spin fluctuations at Q (Stt/S, 2tt/3) due to qlD 
bands are obtained in the three-orbital model. The wavenum¬ 
ber of the incommensurate peak of x® (?) is consistent with 
the neutron measurements. 

We also study the superconducting susceptibility, and the 
obtained TSC gap functions on FSa and FS/3 are shown in 
Fig. ID (c). These obtained results are very similar to those 
in the two-orbital model. The TSC gap function on FS 7 is 
very small since the orbital fluctuations do not develop in dxy- 
orbital. The temperature dependences of the TSC (SSC) sus¬ 
ceptibilities, Xt(^)’ in Fig- E] (d), together with X^{Q) and 
X®(Q). Therefore, the spin and quadrupole susceptibilities 
as well as the TSC gap function obtained in the three-orbital 
Hubbard model are very similar to those given in the two- 
orbital model. 

III. SELF-CONSISTENT VERTEX-CORRECTION (SC-VC) 
METHOD: ANALYSIS OF THE THREE-ORBITAL MODEL 

In this section, we analyze the three-orbital Hubbard model 
by using the SC-VC method with the self-energy correction 
(SC-VCe method) in order to conhrm the validity of 

the results of the RG-tcRPA method shown in Sec. Ill B I In this 
section, we set the charge Stoner factor ^ 0.98 in the SC-VCs 
method. In Fig. | 6 ] we show the obtained (a) spin susceptibil¬ 
ity x^{q)’ (b) quadrupole susceptibility x^{q)’ and (c) mass- 
enhancement factor for FS-(5 (6 = a,/ 3 , 7 ), Zs{k) = 1 — 
Ke-^'Ssik, e) |^_g, in the case of U — 3.05 and J/U = 0.12 
at T = 0.02. We obtain the relation Za^p/Z-^ ~ 0.7, which is 




- 0.02 0 0.02 


FIG. 6: (color online) (a) X*(<7)> (h) and (c) Zs (6 = a, j3, 7) 

given by the SC-VCe method for the three-orbital model without 
SOI. The obtained relation Za^pjZ^ = 0.7 is consistent with the 
result given by the FLEX approximation jsH]. (d),(e) The TSC gap 
functions on a-band and /3-band obtained by the SC-VCe method. 


consistent with the dHvA and ARPES measurements M as 
well as the result of the FLEX approximation m. (The mass- 
enhancement is mainly caused by spin-fluctuations, because 
of the factor 3 in T.{k) ~ J2q + q) id) + 

in the SC-VCs method.) Since the spin fluctuations in the 
7 -band is suppressed by the large (fc), we obtain the exper¬ 
imental strong spin fluctuations at Q Ri (27r/3, 27r/3) due to 
qlD bands 1 I 31 I] . 

A. Gap Equation without the SOI: Verification of the 
orhital-rspin fluctuation mediated TSC nsing the 
SC-VCe method 

Here, we analyze the linearized gap equation. For the triplet 
(a = t) and singlet (a = s) superconductivity, the gap equa¬ 
tion in the absence of the SOI is given as 

q;,/3,7 

AfA^(q) = 

q' 

x|G^,(9')PA^'(g'), (9) 

where Gy {q) is the Green function for band p with self-energy 
correction, and q = {q, e„), where e„ is the fermion Matsub- 
ara frequency. Also, 

vr\q,<^) = X UMUy-l'{-q)Uy,,Mr 

1,1' ,m,m' 

^Uy,,m'{-q'rVt^-^^'^^'{q-q'), ( 10 ) 

where Uy-i{k) = (p, k\l, k) is the unitary matrix. I, I', m, m' 
represent the d-orbital. Va{q) — baV^{q) + CaV°(q) is the 
pairing interaction with the vertex correction for the SC gap 
(A-VC) in the orbital basis, given in Ref. For U = 3.05 
and J/U = 0.12 at T = 0.05, the eigenvalue of the TSC state, 
Af' = 0.304, is larger than that of the SSC state, Af' = 0.291. 
Figures| 6 ](d) and (e) show the obtained TSC gap functions on 
a-band and /3-band at the lowest Matsubara frequency. The 
relation A“’^ ex sinSfcx is satisfied approximately. Without 
the SOI, the 2D band does not contribute to the TSC since the 
relation |A2(/c)|/|A“’^(fc)| ^ 0.1 holds. 

We also performed the SC-VCs study for U = 3.2 and 
J/U = 0.15, and verified that the obtained numerical results 
are similar to Figs.| 6 ](a)-(e) for U = 3.05 and J/U — 0.12. 
Especially, the relationship Af > is satisfied, and the 
(sin3fca;)-like TSC gap function is obtained. 

Thus, we obtained the strong orbital-Hspin fluctuations at 
Q Ri (27r/3,27r/3) as well as the TSC state in the qlD 
bands, by analyzing the three-orbital model using the SC- 
VCs method. These results are consistent with the results 
given by the RG-tcRPA method in Sec.lU 

B. Gap Equation with the SOI: SOI-induced large gap on the 
FS7 due to orhital-mixture 

In this subsection, we discuss the important role of the SOI 
in the superconducting state. If the SOI is neglected, the rela¬ 
tion |A2(fc)|/|A“’^(fc)| < 0.1 holds in both the SC-VCs 
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FIG. 7; (color online) (a) dxy and d^z+yz orbital weights on the 7 -FS 
given by the WIEN2k band calculation with the SOI. Due to the SOI, 
dxz+yz orbitals are mixed in the FS 7 around 6 = it/A. (b) Orbital 
weights on the 7 -FS and (c) three FSs in the three-orbital model with 
SOI (A = 0.4). (d)-(f) The Eu TSC gap functions (fe) oc sin 
(S = a, P, 7 ) obtained by the SC-VCe method. We put A = 0.4 in 
the gap equation. Ay{kz:, ky) = A^{ky, k^P) are not shown, (g) The 
magnitude of the chiral TSC gap |A'^(fe)l = A^(k)'^ + A^{k)^. 


and RGh-cRPA methods. However, the experimental rela¬ 
tion C/T ^ T below Tc {C is the specihc heat) indicates 
that max{|A2(fc)|} takes large value. We expect that large 
A2!(fc) is induced from the gap on the FS/3, due to the “mix¬ 
ing of the orbital character” caused by the SOI. Figure [T] (a) 
shows the orbital weights 

N^y{0)/N2{0) on the FS7 obtained by the WIEN2k with 
SOI, where (0) is the DOS of the Z-orbital (I = xz, yz, xy) 
and iVj(O) = N2z+yz{^) + ^2y{^)- Thus, the SOI-induced 
orbital-mixture is very large, as reported by the hrst-principles 
study and observed by Spin-ARPES measurements 1^ . 
This result means that the induced gap on the 7-band, which 
is approximately given as |A2(fe)| ~ rxz+yz{k)\A2.{k)\, be¬ 
comes comparable to other gaps around 6 ^ njA. 

Here, we introduce the SOI iTsoi = A h ■ •St to 
the present three-orbital model, and put A = 0.4. Eig- 
ure I2] (b) shows the orbital weights N2y{0)/N2{0) and 
^xz+yzi^)/^di^) T'S7 in the present three-orbital 

model with A = 0.4. The ESs for A = 0.4 are shown in 
Fig- 13 (c). By taking the SOI into account, the linearized gap 
equation in terms of the SC-VCs method for a = t, s is given 


as 


EC;.';P3P4(9.'?') 

q' m' P3,p4 

x|G^,(9')PACaP4(9'), (11) 

where y, y! represent the bands with SOI, and pi is 
the pseudo-spin (fl", JJ.) that represents the Kramers doublet. 
Aa P1P2 (9) 1^ even (odd) with respect to the exchange of the 
pseudospin for a = t (s). Using the d-vector, the TSC gap 
functions are expressed as 


(9) 

(9) 

f +id>2{q) d^iq) 

{ d2{q) d2{q)+id>2{q) 


At,ti4(9) / 


A'" 


. ( 12 ) 


The pairing interaction is given as 


PlP2',P3P4 


(9,9') 


U„, 


^j '-^ppi;i(Ti 

1 , 1 ' ,m^m' ,<7i 

-.'A* 


{q)U, 


PP2;i'o-2 


(- 9 ) 


xU, 


xK 


(9-9'): 


/T/03;mf73 vy ) ^ fip4‘,rn'a4 
l,m;m' , 1 ' 

<7i(T3;cT4<72 


i-Q'r 


(13) 


where Uyp-^ia{k) = {yp,k\la,k) is the unitary matrix {I 
is the d-orbital, cr is the real spin). We approximate that 
lAiCT2;o-3cr4(9) is given by the SC-VCs method without the 
SOI, introduced in Ref. lE^I . (Note that 14 i(T 2 ;ct 3 0-4(?) = 

-iy‘^(,q) 5 a,,aJai,a 3 +V’^{q)aa,,a 2 ' ^^4 .<t 3)/2, where 
is the interaction for the charge (spin) channel with A-VC 
||23] .) This approximation is reasonable since the pairing in¬ 
teraction is induced by higher-energy processes, whereas the 
SOI is not important except near the band crossing points. In 
fact, the relation xHq) ~ X%[y){Q.) is obtained in the RPA 
for A = 0.4 even for the spin Stoner factor as = 0.99. 
Thus, spin-fluctuation-driven VC will give strong orbital fluc¬ 
tuations even in the presence of the SOI. 

The six-fold degeneracy of the TSC state is lifted by the 
SOI. Eor U = 3.05 and J/U = 0.12 at T = 0.05, the two- 
dimensional Eu TSC state (dx = dy = 0, dz = ^x, Ay) 
is obtained with the largest eigenvalue A® = 0.237. (A^ = 
0.234 for the Aig SSC state, and A^ = 0.232 for the helical 
Aiu TSC state; dx = A^,, dy — Ay, dz = 0.) The Eu TSC 
gap functions at the lowest Matsubara frequency are shown in 
Eigs.|3(d)-(f) in the presence of the SOI with A = 0.4. 

In the case of TSC state, the chiral superconductiv¬ 
ity dz = Ax + iAy is expected to be realized below T^. 
Eigure 0 (g) shows the magnitude of the chiral gap state 

|A'^(fe)| = A2(kY + Ay(kY for J = a,(3,x- There¬ 
fore, relatively large gap on the ES7 is induced by that on 
the ES/3 around 0 ~ 7r/4 due to the SOI, which is approxi¬ 
mately given as |A'*'(A;)| ^ rxz+yz{k)\A^(k)\. This result is 
consistent with the experimental report by the field-orientation 
dependent specific heat measurement . It is our important 
future problem to explain the approximate T-linear specific 
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heat. Also, we will perform the RG+cRPA analysis for the 
three orbital Hubbard model with SOI in future. 

In contrast, if the TSC is mainly realized on the FS7, the 
gap on the FSa induced by the SOI should be very small, 
since the orbital mixing is very small in the a-band. In this 
case, the “residual specific heat” appears for T ^ Tc due to 
the tiny gap on the FSa, although it is inconsistent with the 
specific heat measurement. 

To summarize this subsection, it is found that substantial 
superconducting gap on FS7 is induced from the FS/3 due 
to the large SOI of Ru-atom (proximity effect), since SOT 
induced orbital-mixture is very large between FS/3 and FS7. 
We also briefly discussed the d-vector ifT^l^ls^ in the pres¬ 
ence of the SOI, which is closely related to the important topo¬ 
logical properties of the TSC state ll37l - l^ . 

IV. SUMMARY 

In the present paper, we proposed that the orbital-tspin 
fluctuation-mediated TSC is realized in Sr2Ru04 by applying 
the RGh-cRPA method and the SC-VCe method to the multi¬ 
orbital Hubbard models. Thanks to the VC that is neglected 
in the RPA, strong orbital and spin fluctuations si q ^ Q 
emerge in the qlD FSs. Then, the TSC is easily induced by 
the coexisting orbital and spin fluctuations. The present the¬ 
ory naturally explains why the TSC overwhelms the SSC in 
Sr2Ru04 despite the strong AFM fluctuations. 


We also analyzed the effect of the SOI on the superconduct¬ 
ing gap function. It is found that the large gap on the FS7 is 
induced from that on the FS/3 due to the SOTinduced “orbital- 
mixture” between FS/3 and FS7. This effect will be important 
to understand the experimental approximate T-linear behav¬ 
ior of C/r below Tc. If the TSC were mainly realized on the 
FS7, the “residual specific heat” should appear for T <C Tc 
since the SOTinduced orbital-mixture between FSa and FS7 
is very small. It is our important future problem to obtain the 
phase diagram of the d-vector in the TSC state. 

Note added in proof. Recently, we became aware of the 
resonant inelastic X-ray scattering study lEH, in which the 
SOTinduced large orbital mixture in Sr2Ru04 discussed in 
the present paper had been confirmed experimentally. 
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